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The holographic dual of a finite-temperature gauge theory with a small number of flavours typi- 
cally contains D-brane probes in a black hole background. At low temperature the branes sit outside 
the black hole and the meson spectrum is discrete and possesses a mass gap. As the temperature 
increases the branes approach a critical solution. Eventually they fall into the horizon and a phase 
transition occurs. In the new phase the meson spectrum is continuous and gapless. At large iV c and 
large 't Hooft coupling, this phase transition is always of first order, and in confining theories with 
heavy quarks it occurs at a temperature higher than the deconfinement temperature for the glue. 

PACS numbers: 



Introduction: The gauge/gravity correspondence is a 
powerful tool to study non-perturbative physics of gauge 
theories in diverse dimensions. The classical supergravity 
regime corresponds to the large- N a , strong 't Hooft cou- 
pling limit of the gauge theory. This allows the study of 
a large class of theories that share some of the important 
features of four-dimensional QCD, such as confinement, 
chiral symmetry breaking, thermal phase transitions, etc. 
In principle, because of its asymptotic freedom, QCD it- 
self is not in this class. This means that calculations of 
certain quantitative properties of QCD, such as the de- 
tailed mass spectrum, for example, will require going be- 
yond the supergravity approximation. However, this does 
not exclude the possibility that some aspects of QCD 
can be studied in this approximation: Some predictions 
of the gauge/gravity correspondence may be universal 
enough as to apply to QCD, at least in certain regimes. 
A suggestive recent example is the gauge/gravity calcula- 
tion of the shear viscosity in the hydrodynamic regime of 
strongly coupled finite-temperature gauge theories. The 
viscosity /entropy ratio is universal for a large class of 
gauge theories in the regime described by their gravity 
duals [1| which, for high enough a temperature, generi- 
cally contain a black hole horizon. Moreover, this ratio 
appears to be surprisingly close to that inferred from ex- 
periments at the Relativistic Heavy Ion Collider (RHIC) 
|2|. It is therefore important to establish as many uni- 
versal features of the gauge/gravity correspondence as 
possible. 

The quarks in QCD transform in the fundamental rep- 
resentation. For a large class of gauge theories, a small 
number of flavours of fundamental matter, N { N c , may 
be described by probe D-branes in the appropriate gravi- 
tational background 3] . At sufficiently high temperature 
T, this contains a black hole 4j. The purpose of this pa- 
per is to exhibit some universal features of this system 
that only depend on these two facts. In particular, we 
demonstrate the existence of a first order phase transition 
for the fundamental matter, as follows: 

At sufficiently small T/Af q (where M q is the quark 
mass) , the brane tension is sufficient to overcome the at- 
traction of the black hole and hence the branes lie outside 



the horizon in a 'Minkowski' embedding — see fig. ^ and 
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FIG. 1: Profiles of D7-brane embeddings in a D3-brane back- 
ground. The thick black circle is the horizon (p = 1). 

below. In this phase the meson spectrum (i.e., the spec- 
trum of quark-antiquark bound states) is discrete and 
possesses a mass gap. At sufficiently large T/Af q , the 
gravitational attraction overcomes the brane tension and 
the branes fall into the horizon yielding a 'black hole' em- 
bedding. In this case, the meson spectrum is continuous 
and gapless. In between, a limiting, critical solution ex- 
ists. We will show that the phase diagram in the vicinity 
of this solution exhibits a self-similar structure. While 
this structure went unnoticed, the phase transition that 
occurs as T/M q increases from small to large values was 
observed in two specific models H, E|- In fact, as was 
first noted in Q for a D6-brane in a thermal D4 back- 
ground, the transition is of first order (in the approxima- 
tions stated above). Rather than dropping continuously 
through the critical solution, the probe brane jumps dis- 
continuously from a Minkowski to a black hole embed- 
ding at some T — T tund . This leads to a discontinuity in 
several field theory quantities, such as, for example, the 
quark condensate (f/'V') or the entropy density. In the 
following, we will see that the critical behavior and, as a 
result, the first order transition are essentially universal 
to all Dp/Dq systems. 
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The case of confining gauge theories is particularly in- 
teresting because, for sufficiently heavy quarks, two dis- 
tinct phase transitions occur . At T = T dacaat the grav- 
itational background changes from a horizon-free back- 
ground to a black hole background. At this temperature 
the gluons and the adjoint matter become deconfined Q, 
but the branes remain outside the horizon and hence sta- 
ble quark-antiquark bound states still exist in a range 
Tdccont < T < T tuad . At T — T fund the phase transition 
described above for the fundamental matter occurs. 
The Setup: The black Dp-brane metric in the string- 
frame takes the form 



ds 2 



H~2 (-fdt 2 + dxi 



du 2 

T 



u 2 dn 2 8 _ p 



(i) 

where H{u) = {L/u) 7 - p , /(«) = 1 - (u a /u) 7 - p and L 
is a length scale (the AdS radius in the case p = 3). 
There is also a non-trivial dilaton = gsH^P" 4, and 
a Ramond-Ramond field Coi... p = if -1 . The horizon lies 
at u = uq. As usual, regularity of the Euclidean section, 
obtained through t — ► it E , requires that t E be identified 
with period 
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T 



47rL 
~T~p 



(2) 



According to the gauge/gravity correspondence, string 
theory on the background above is dual to a (p + 1)- 
dimensional supersymmetric gauge theory at tempera- 
ture T. In some cases one periodically identifies some of 
the 'Poincare' directions x v in order to render the the- 
ory effectively lower-dimensional at low energies; a proto- 
typical example is that of a D4-brane with one compact 
space direction. Under these circumstances a different 
background with no black hole may describe the low- 
temperature physics, and a phase transition may occur 
as T increases Q . In the gauge theory this is typically a 
confinement /deconfinement phase transition for the glu- 
onic (or adjoint) degrees of freedom. Throughout this 
paper we assume that T is high enough, in which case 
the appropriate gravitational background is always l|T]l. 

Consider now a Dq-brane probe that shares d spacelike 
'Poincare' directions with the background Dp-branes 
and wraps an S n inside the 5 8 ~ p . We will assume that 
the Dq-brane also extends along the radial direction, so 
that q = d + n + 1. In the gauge theory this corresponds 
to introducing fundamental matter that propagates 
along a (d + l)-dimensional defect. To ensure stability, 
we will assume that the Dp/Dq intersection under 
consideration is supersymmetric at zero temperature. 
Under these conditions the RR field sourced by the 
Dp-branes does not couple to the Dq-branes. Two cases 
of special interest here are the D3/D7 (n — 3) 6] and the 
D4/D6 (n = 2) @ systems. If one of the D4 directions 
is compact, then both cases can effectively be thought of 
as describing the dynamics of a four-dimensional gauge 
theory with fundamental matter. 



Criticality and Scaling: In this section we fol- 
low Q closely. We begin by studying the behaviour of 
the brane probe near the horizon. In order to do so we 
write 



dtt 2 8 _ p = d9 2 + sin 2 9 dill + cos 2 6 dn 2 7 _ p _ n 



(3) 



u = uq + ttTz 



p-3 

L \u J 



with T the temperature above. Expanding the metric 
to lowest order in z, y gives Rindler space together with 
some spectator directions: 



ds 1 



-(2TrT) 2 z 2 dt 2 +dz 2 +dy 2 +y 2 dn 2 n 



-dx 2 y 



(4) 



The Dq-brane lies at constant values of the omitted co- 
ordinates, so these play no role in the following. The 
horizon is of course at z = 0. The Dq-brane embedding 
is specified by a curve (z(a),y(a)) in the (z, y)-plane. 
Since the dilaton approaches a constant near the hori- 
zon, up to an overall constant the Dq-brane (Euclidean) 
action is simply the volume of the brane, namely 



/ cx J da\/I 2 ~ 



y 2 zy r 



(5) 



where the dot denotes differentiation with respect to a. 
This is precisely the action considered in ref. Q- In the 
gauge z — a the equation of motion takes the form 



z yy + (yy - nz ){ 1 + y 2 ) = o • 



(6) 



Solutions fall into two classes that we call 'black hole' em- 
beddings and 'Minkowski' embeddings; see fig. ^ Black 
hole embeddings are those for which the brane falls into 
the horizon, and may be characterised by the size of the 
induced horizon, y . The appropriate boundary condi- 
tion in this case is y = 0, y = yo at z — 0. Minkowski em- 
beddings are those for which the brane closes off smoothly 
above the horizon. These are characterised by the dis- 
tance to the horizon, zq, and satisfy the boundary condi- 
tion z = 0, z = zq at y = 0. The limiting solution is the 
critical solution y = y/nz, which touches the horizon at 
the point y = z = 0. 

The equation of motion @ enjoys a scaling symmetry: 
If y = f(z) is a solution, then so is y — f(pz)//j for 
any real positive fi. This transformation rescales z$ — > 
Zq/h for Minkowski embeddings, or yo yo/ (J- for black 
hole embeddings, which implies that all solutions of a 
given type can be generated from any other one by this 
transformation. 

Consider now a solution very close to the critical one, 
y(z) — \fnz + £(z). Linearising the equation of motion 
©, one finds that for large z the solutions are of the 
form = z v ±, with u± = -n/2 ± ^n 2 - 4(n + l)/2. 
If n > 5 these exponents are real, whereas if n < 4 they 
have non-vanishing imaginary parts that lead to oscil- 
latory behaviour. Under the assumptions stated above 
n > 5 implies d < 1 (an example of this is a D3/D7 in- 
tersection over a string). Since we are mostly interested 
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in higher dimensional defects, we will henceforth assume 
that n < 4. In this case it is convenient to write the 
general solution as 



'nz+ 



T- 



asin(alogTz)+6cos(o;logT2;) , (7) 



where a = y 4(n + 1) — n 2 /2 and a, b are dimensionless 
constants determined by z$ or j/o- It is easy to show 
that under the rescaling discussed above, these constants 
transform as 



1 



cos(alog/x) sin(alog/i) 
- sin(a log p) cos(a log p) 



(8) 



This transformation law implies that the solutions ex- 
hibit discrete self-similarity and can be used to derive 
critical exponents that characterise the near-critical be- 
haviour. We refer the reader to 0, @ for details but we 
emphasise that this behaviour depends only on the di- 
mension of the sphere, and is therefore universal for all 
Dp/Dq systems with n < 4. 

Each near-horizon solution gives rise to a global 
solution when extended over the full spacetime. Each of 
these solutions is characterised by a quark mass M q and 
a quark condensate (ipip) that can be read off from its 
asymptotic behaviour. Both of these quantities are fixed 
by zq or y . As we will see, the values corresponding 
to the critical solution, M* and (^ip)*, give a rough 
estimate of the point at which the phase transition 



Fundamental Phase Transitions: In order to 
study the global solutions, it is convenient to introduce 
an isotropic, dimensionless radial coordinate p through 



(uop) 



^£ 7 -P I i 7-r> 



(9) 



Note that the horizon is at p = 1. Just for concrete- 
ness, we will now assume that the Dp/Dq system under 
consideration is T-dual to the D3/D7 one, in which case 
(p — d) + (n + 1) = 4. Under these circumstances, the 
Euclidean Dq-branc action density in the black Dp-brane 
background is 



AT 



/"*(-) 

J Putin \ U 0Pj 



d-3 



1 



1 



5 2(7-p) 



x(l-X 2 ) 2 y/1 - X 2 + P 2 X 2 



where Af is a normalisation constant: 



Af = 



AT 



(10) 



(11) 



T Dq = l/(2ir£ s ) q g s £ s is the Dq-brane tension, f2„ is the 
volume of a unit n-sphere, x = cos 9, and x = d-x/dp. 
Multiple flavours arise by introducing N t coincident 
probe branes. Up to a numerical constant the normali- 
sation factor is found to be 



(12) 



where g 2 lf (T) = \T P ~ 3 is the dimensionless effective 't 
Hooft coupling for a (p + l)-dimensional theory at tem- 
perature T 0, A = g 2 M N c and we have used the stan- 
dard gauge/gravity relations g 2 M ~ g s ^~ 3 and L 7 ~ p ~ 
gsN^il'P with l s the fundamental string length. 

The equation of motion that follows from (|10fl leads to 
the large-p behaviour [21) 



m c 

X = — + ~ 

p p n 



(13) 



The constants m, c are_related to the quark mass and 
condensate through 



M n = 



worn 

9-p 

2— tt£? 



2ire 2 s (n - i)n n r DqM gc 



(14) 

This implies the relation m^ 5 p )' 2 = M /T between the 
dimensionless quantity m, the temperature T and the 
mass scale 



M = 



7- 



P 



(15) 



This is precisely the scale of the mass gap in the discrete 
meson spectrum at temperatures well below the phase 
transition |{J fToL ITU . We see here that it is also the 
scale of the transition temperature for the fundamental 
degrees of freedom, T fund ~ M, since this takes place at 

TO ~ 1. 

The key observation [S| is that the values (to, c) of 
a near-critical solution are linearly related to the cor- 
responding values in the near-horizon region. Combin- 
ing this with the transformation rule (jHJ) for the near- 
horizon constants (a, 6) and eliminating p, we deduce 

that (to — to*)/z 2+1 and (c— c*)/z 2+1 are periodic func- 
tions of (a/2n) log zq with unit period for Minkowski em- 
beddings, and similarly with zq replaced by yo for black 
hole embeddings. This is confirmed by our numerical 
results, as illustrated in fig. [21 

The oscillatory behaviour of to and c as functions of zq 
or yo implies that the quark condensate is not a single- 
valued function of the quark mass. The preferred solution 
will of course be the one that minimises the free energy 
density of the Dq-brane, F — T/ Dq . This quantity con- 
tains a volume divergence, as can be seen by using the 
asymptotic behaviour l|13l) . It therefore needs to be reg- 
ularised and renormalised. We achieve the former by 
replacing the upper limit of integration by a finite ultra- 
violet cut-off p max . The latter can be done by subtracting 
the free energy of a fiducial embedding, as was done for 
the D4/D6 system in ref. 0. The asymptotic geom- 
etry of a D7-brane in a D3 background is AdSs x S 3 , 
so in this case the more elegant method of holographic 
renormalisation f° r brane probes is available. 
This consists of adding to the brane action the boundary 
'counter-term' 



Af 



2\2 



Amcj . 



(16) 
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FIG. 2: Quark mass as a function of the distance to the hori- 
zon zo for D7-brane Minkowski embeddings in a D3-brane 
background. 

The renormalised brane energy / = / Dq + I hound is then 
finite as the cut-off is removed, p max — > oo. The results 
for the D3/D7 case are shown in fig. We see that 
the transition occurs discontinuously by jumping from a 
Minkowski embedding (point A) to a black hole embed- 
ding (point B). We emphasize again that this first order 
transition is a direct consequence of the multi- valued na- 
ture of the physical quantities brought on by the critical 
behaviour described in the previous section. It may be 
possible to access this self-similar region by super-cooling 
the system. 

It is interesting to ask if the strong coupling results 
obtained here could in principle be compared with a weak 
coupling calculation. It follows from our analysis that the 
free energy density takes the form F = JVT f (m 2 ) , where 
the function / can only depend on even powers of m 
because of the reflection symmetry x ~ > ~ X- The strong 
coupling limit g ot{ — ► oo corresponds to m — > 0, which 
may be equivalently regarded as a zero quark mass limit 
or as a high-temperature limit. In this limit the brane 
lies on the equatorial embedding \ — an <i slices the 
horizon in two equal parts. In general /(0) is a non-zero 
numerical constant; in the D3/D7 case, for example, a 
straightforward calculation yields /(0) = —1/2. Thus at 
strong coupling the free energy density scales as 

F ~ N t N T d+1 g etf (T) ^ . (17) 

The temperature dependence is that expected on di- 
mensional grounds for a <i-dimcnsional defect, and the 
N f N c dependence follows from large- TV counting rules. 
However, the dependence on the effective 't Hooft cou- 
pling indicates that this contribution comes as a strong 
coupling effect, without direct comparison to any weak 



coupling result. The same is true for other thermody- 
namic quantities such as, for example, the entropy den- 
sity S = —dF/dT. We remind the reader that the back- 
ground geometry makes the leading contribution to the 
free energy density [Tjj 

F ~ N^T^g^T)^ 1 , (18) 

which corresponds to that coming from the gluons and 
adjoint matter. Of course, for p = 3, the effective cou- 
pling factor is absent and this bulk contribution differs 
from the weak coupling result by only a factor of 3/4 [l5| . 

Discussion: We have shown that, in a large class 
of gauge theories with fundamental matter, quark- 
antiquark bound states survive the deconfining phase 
transition for the gluonic degrees of freedom provided 
M > T doconf , where M is the typical mesonic scale. 
This is potentially interesting in connection with QCD, 
since in QCD heavy quark mesonic bound states with 
M ^> T doconf ~ 175 MeV certainly exist. One generic 
feature of the low-lying mesons in the class of theories 
discussed here is that they are extremely deeply bound 
at strong coupling as is apparent from eq. JTSJl. 

It is intriguing that the mesonic states claimed to ex- 
plain certain features (such as the entropy density) of the 
strongly coupled quark-gluon plasma formed at RHIC are 
also deeply bound fl6| . It would be remarkable if a pre- 
cise relationship between the two could be established. 

Thermal phase transitions in gauge theories with spon- 
taneously broken chiral symmetries at zero temperature 
are particularly interesting, since they raise the ques- 
tion of whether chiral symmetry is restored at the phase 
transition. This has been found to be the case [l7| in 
the D4/D8/D8 holographic model of 18]. In this model 
the D8/D8 pair is connected before it falls into the hori- 
zon, but splits into two disconnected pieces when it does. 
At this point the symmetry is enhanced from the diago- 
nal U(Nf) v subgroup to the full £7(iV f ) L x U(N f ) R group. 
This contrasts with the cases discussed in this paper, in 
which the branes remain connected when they fall into 
the horizon. From the gauge theory viewpoint, this dif- 
ference in topologies is due to the fact that in the first 
case the fundamental matter lives on a defect /anti-defect 
pair, whereas in the second it lives on a single defect. 

The detailed results found in this paper rely on the ap- 
proximation that 1/A C , 1/A — ► with AT f fixed. However, 
the fact that the phase transition is first order implies 
that it should be stable under small perturbations, and so 
its order and other qualitative details should hold within 
a finite radius of the 1/A C ,1/A expansions. Of course, 
finite- N c and finite-A corrections may eventually modify 
the behaviour uncovered here. For example, at large but 
finite N c the black hole will Hawking-radiate and each bit 
of the brane probe will experience a thermal bath at a 
temperature determined by the local acceleration. This 
effect becomes more and more important as the lower 
part of a Minkowski brane approaches the horizon, and 
may potentially blur the self-similar, scaling behaviour 
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FIG. 3: Quark condensate and free-energy density for a D7 in a D3 background; note that JV oc T 3 . The blue dashed 
(red continuous) curves correspond to the Minkowski (black hole) embeddings. The dotted vertical line indicates the precise 
temperature of the phase transition. 



found here. Note that this effect is of order l/iV c 2 , and 
therefore subleading with respect to the order- N f /N c cor- 
rection to physical quantities from the presence of the 
brane probes. 

Finite 't Hooft coupling corrections correspond to 
higher-derivative corrections both to the supergravity 
action and the D-brane action. These may also blur 
the structure discussed above. For example, higher- 
derivative corrections to the D-brane equation of motion 
are likely to spoil the scaling symmetry of eq. © , and 
hence the self-similar behaviour. These corrections also 
become important as the lower part of a Minkowski brane 
approaches the horizon, since the (intrinsic) curvature of 
the brane becomes large there. 

Yet another type of correction one may consider is due 
to the backreaction on the background spacetime of the 
Dq-branes, whose magnitude is controlled by the ratio 
N t /N c . These have been considered for the D2/D6 sys- 
tem in ref. 19' . In particular, this ref. finds that the en- 
ergy density scales as F ~ N^ 2 N^ 2 T 3 , which obviously 



differs from (|17|l with p = 2, d = 2. This discrepancy 
is not at all a contradiction, and has the same origin as 
the discrepancy found for the meso n sp ectrum [lfj . This 
is the fact that the calculation in [H| applies in the far 
infrared of the gauge theory, whereas that presented here 
applies at high temperatures, i.e., at T > g%. . 
We will return to these and other issues in |20| . 
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